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I. INTRODUCTION* 1 ' 2 ) 


The study of point defects in semiconductors has a long and honorable history (see 
Table I). In particular, the detailed understanding of shallow defects in common 
semiconductors traces back to the classic work of Kohn and Luttinger. However, the 
study of defects in narrow gap semiconductors represents a much less clear story. 
Here, both shallow defects (caused by long range potentials) and deep defects (from 
short range potentials) are far from being completely understood. 

The study of narrow gap semiconductors (see Table II) is important to NASA 
because typical ones such as Mercury Cadmium Telluride and Mercury Zinc 
Telluride are difficult to grow uniformly under the stress of gravity induced 
convection. Also, the narrow gap semiconductors are useful for infrared detection 
which in turn is useful for remote sensing and a variety of military purposes. 

One aspect of studying crystal perfection is to study the natural and artificial defects 
which appear in these materials. In our study we calculate the properties of 
numerous examples of deep defects and discuss some problems associated with the 
study of shallow defects (see Table HE). All defects we study are point defects. Until 
one can identify and characterize defects, it is impossible to make clear statements 
about crystal perfection. 

A full characterization of crystal perfection involves not only a theoretical study of 
all common defects but also experimental verification of the models. The simplest 
model of shallow defects assumes they are fully ionized. There is not a lot one can 
do with such a model. However, refinements tend to get one into many 
complications, as we discuss in the appendix. Donor defects can perhaps be best 
studied by photoluminescence experiments. It is also of interest that the scanning 
tunneling microscope, which can be used to measure the local density of states, can 
also be used to study deep defect energy levels. 

In this study, all results are calculational and our focus is on the chemical trends of 
deep levels in narrow gap semiconductors (see Table IV). We study substitutional 
(including antisite), interstitial and ideal vacancy defects. For substitutional and 
interstitial impurities, the effects of relaxation are included. For materials like 
Hgj^Cc^Te, we study how the deep levels vary with x, of particular interest is what 

substitutional and interstitial atoms yield energy levels in the gap i.e. actually 
produce deep ionized levels (see Table IV). Also, since the main technique utilized 
is Green's functions, we include some summary of that method. 
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n. DEEP LEVELS OF POINT DEFECTS IN SEMICONDUCTORS 


A. Green's Function Method^' 6 ) 

Let us suppose that an atom of the perfect crystal is replaced by a defect. If the 
Hamiltonian of the perfect crystal is H Q and if the change in potential produced by 

the defect is V, then Schrodinger's equation becomes 

(H 0+ V)|>P>=e|4'>. (!) 

The host crystal's Green's function is defined as 

G°|E) = (E-H 0 )‘ 1 . (2 ) 


Combining these two equations, we have 
(i - G° (E) V ) IT>=0. 


(3) 


The defect energy levels E in the band gap of a semiconductor are thus given by non- 
trivial solutions of the determinantal equation: 

Det (i - G° (E) V ) = 0 . ( 4 ) 

For a point defect in a tetrahedral site of a zinc blende material, the point group is 
T d . Such a defect can have both non-degenerate Aj (s - like) state or triply 
degenerate T 2 (p - like) states. Assuming these states, Eqn (4) can be factored as 



_ _ 

3 

r l 

det 

I - G°(E) V 

= n det 

Ii-G-Vi 


i =0 


(5) 


for i = 0 corresponding to A| states and i = 1, 2, 3 to T 2 states. 
For Aj states only 


det 


i s <v s 


= 0 


( 6 ) 
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and for the T 2 states 


det 


V g p v pj 


= 0 , 


(7) 


for suitably defined G s , V s and G , V 


U s = 1 /s andU D = 1 /p 

F 'Jaa (see Appendix B4), the deep levels in 


Using the equations ^aa 

the band gap can be determined as a function of U g and Up , which is what we mean 

by impurity potential in Figs. 3-6. Here U s = P s (E kl ^ > - E hos J an d 


Up - P p (E^p- E^ ost J t where E irnp and E host are the ground state orbital energies of 
defect and host atom, respectively. The relationship between U g , Up and V is 


explained in the appendix. Since the Green's function G°(E) is a function of energy 
E, we can calculate G°(E) for E in the band gap. Thus we can obtain a plot of U g vs E 
or E vs U g (and similarly for Up). For each substitutional defect, the difference in 

energies will be determined as will P s , (Pp). Thus we can locate the deep defect 
energy levels. 
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B. Comments on Ideal Vacancy Model( 7 ) 

Vacancies are the most common native defects. The ideal vacancy model is an often 
used model for predicting deep vacancy levels. This model is described in detail in 
the paper by Das Sarma and Madhukar. The model assumes that a vacancy is 
formed by removing an atom from the crystal, but leaving all other atoms in the 
same position. Das Sarma and Madhukar proved that this is equivalent to setting 
U s and Up equal to infinity. From our previous equations, this is equivalent to 

setting 


GL(E) = 0 


gL(E) = o 


This, in turn, is equivalent to setting atomic orbital energies equal to 


(8) 
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C. Band Structure Theory of II- VI Semiconductors^ 8 ) 

Unlike III-V semiconductors, II- VI semiconductors may have strong spin-orbital 
coupling, which must be included in the band structure calculations. Based on Vogl 
et al. band structure theory, Kobayashi proposed a theory for II- VI semiconductors 
which included spin-orbital coupling. The Hamiltonian for the perfect crystal in 
this theory is (assume that the defect is on the anion site) 


H 0 = X liaao>Ej<iaaol +1 icad>Ej<icadl 


+ X 


liaao>V ii <jc adl + h.c. 


+ H sn , 


( 9 ) 


where a is the spin, o specifies the unit cell, and d is the nearest-neighbor vector, 
is the spin-orbit Hamiltonian, which can be written 


H S „ = X 


I i a cto >2X a L a 


a a < j aaol+ I icad>2A, c L c 


a c < j cadi 


J (10) 


where and are the spin-orbit parameters. As usual, the band structures of the 
alloys were obtained within the virtual crystal approximation. 

However, for most point defects that we are interested in, the spin-orbit interaction 
is weak and can be ignored. Then, including the degeneracy due to spin just doubles 
the degeneracy of each level. This makes the states doubly degenerate, and the T 2 

states six fold degenerate. For this case, the previous Hjalmarson et al theory is still 
valid. 
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D. Lattice Relaxation and Molecular Dynamics Approach^ 9 " 13 ) 

1. When the off-diagonal elements of the defect potential are equal to zero, the 
effects of the lattice relaxation are excluded since only these elements are 
related to defect-bond length. To include these effects, we must include the 
off-diagonal elements of the defect potential. Using the nearest-neighbor 
approximation, the form of the defect potential is 

V =SV i = I i I iao>V i <iaol 

+ I y liao>a ij <jcdl (11) 

where i can be s or p, and <Xjj define the off diagonal elements. When all the 
symmetries are built in, we find the form of the matrix element Vj is 


uL a a a a 

a a 


a 

V i= a O 

a 
a 


The determinantal equation then leads to the following: 



G l,+ G lc , 4GL,G‘ c -G^G;,-3Gl.G g lr 2 

_i i „i * 


( 12 ) 


( 13 ) 


where G a is the elements of Green's function matrix coupling two nearest- 
neighbor host atoms. This gives an implicit equation for the i - symmetric deep 
levels E as a function of Uj and oq. 


The form of a can be inferred from arguments due to Harrison which imply that the 
nearest neighbor matrix elements scale as the inverse bond length squared. We thus 
find 
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( 14 ) 



where dj is the bond length between the defect and one of the nearest neighbors, and 
d H is the bond length of the host crystal, with Cj being a proportionality constant. 

Since otj is related to dj it is also directly related to the lattice relaxation. Further the 
sign of dj determines the direction of lattice relaxation. 

If we can determine a for a particular defect, the deep levels induced by this defect 

can then be predicted by the previous equation. In order to determine a, we need to 
calculate dj and d H . For well known semiconductors d H is known. The 

determination of dj is more difficult, and we now discuss this. 

2. Molecular Dynamics Approach 

(a) In order to calculate a for a particular defect in a host, the defect bond 
length must be specified . A simple model will first be presented, and 
then the more general molecular dynamics approach. The simple 
model uses the concept of covalent radius. It is assumed that 

dj = f i + r H (15) 

where rj and r H are the covalent radii of the defect and the nearest-neighbor atoms 

respectively. Such radii are tabulated for most atoms of interest in the present study. 
The host bond length is given by 

a H -r H +r H' (16) 

For a particular defect in a host, if dj > d H we have outward relaxation. If dj < d H , 
we have inward relaxation. 

This simple model works very well for most group IV and III-V semiconductors. 
However, it is not suitable for II- VI semiconductors because they are not purely 
covalent materials, but a mixture of covalent and ionic bonds. For these we must 
use a molecular dynamics approach. 

(b) For an ideal defect, the four nearest-neighbor host atoms will remain 
in their perfect-crystal positions. However, in a more realistic 
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treatment, these atoms experience a net force and the lattice distorts to 
a new configuration where that force is zero. The total force on one of 
the nearest-neighbor atoms can be divided into repulsive and attractive 
parts, which can be represented as 


F =F r + F a 

X X X 


(17) 


where x is along the defect-host bond direction. The origin of F x is the repulsion 

between electrons in overlapping states, and this force can be computed from the 
empirical relation 



where A is a constant determined by the condition 


Fx(d H ) = F*(d H ) 


(19) 


3V 

which assumes the same functional form for 3x for the host crystal as for the host 

a 

crystal with defect. The attractive force F x can be computed by use of the equation 



( 20 ) 


where £p is the Fermi energy, G is the Green's function matrix including the defect 
and V is the defect potential (see Eqn. B38). 

Using these two equations and beginning with the nearest neighbor host atoms in 
their perfect-crystal positions, we can calculate the motion of each of these atoms for 
a small time interval. At, using Newton's second law with the standard molecular 
dynamics approach. One lets the process continue until the net force approaches 
zero. Thus dj is determined. 
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E. Interstitial Defects - Brief Comments ( 14_15 ) 

Self interstitial impurities and vacancies are the most important native defects in 
semiconductors. In II- VI materials, it is believed that a substitutional group - 1 
acceptor may spontaneously move to an interstitial site. 

For a substitutional defect in zinc blende structure semiconductors, the defect has a 
tetrahedral symmetry in the materials. That is, there are four nearest-neighbor 
atoms and twelve second-nearest-neighbor atoms around the defect. If the distance 
between the defect and the nearest neighbor is d, then the distance between the 
defect and its second neighbor is 1.64 d. The situation of tetrahedral-site interstitial 
impurities is different. In this case, there are four nearest neighbors around an 
interstitial impurity with distance d’, but only six second neighbors around it with a 
slightly greater distance 1.15 d’. Hence, the nearest-neighbor approximation, which 
is often used for substitutional defects, is not acceptable for interstitial impurities. 

We will calculate results for deep levels introduced by interstitial impurities in 
MCT, MZT, and MZS. The theory we will use will take into account the above 
complications but will be somewhat similar to what we have done before. We will 
give our results in the final report. 
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F. Summary of Results for MCT, MZT and MZS for Substitutional Defect 

In Table 4, we summarize the results of several calculations we have done including 
both substitutional vacancy and antisite defects. In the final report, we will include 
these, as well as results in interstitial defects. Here we just give some typical results. 
In Table 5 we give a few typical examples of the calculations including the effects of 
lattice relaxation. 

In Figure 1 (for MZT), we show several possible cation-site s-like (A x ) levels which 
may be formed by (for example) Te (Antisite), C, At and I as a function of 
composition x. The vacancy level below the valence band edge is considered 
because it is possible that lattice relaxation or other effects may raise it into the 
bandgap. In Figure 2 (For MZS), we see anion site, p-like (T 2 ) levels due to Cd, Zn 
(antisite). Mg and Be as function of composition x. 

In Figs. 3, 4, 5 we summarize some results including the effects of lattice relaxation 
for anion-site p-like deep levels which we can get for MZT (x = 0.15), MCT (x = 0.22) 
and MZS (x = .08). In each case the top of the valence band is at E = 0.0 and the 
bottom of the conduction band is at E = 0.1. In Fig. 3, we see that Zn may have an 

impurity level in the band gap for a = 0 (no relaxation). Similar comments may be 
made about Cd in MCT and Be in MZS. In each of these particular cases relaxation 
may take the level out of the band gap or deeper in the band gap (see Table 5). In Fig. 
6, the ideal vacancy levels are modeled by setting the impurity energy = ± «> (or - 120 
is nearly the same). We see the vacancy levels appear below the top of the valence 
band for this case. For other situations we do find some vacancy levels in the band 
gap for MZS. We will provide a detailed summary of our results in the final report. 
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m. SHALLOW IMPURITY LEVELS IN II- VI MATERIALS 
A. General Discussion* 16 ' 17 * 

We give , in the appendix, a fairly complete description of shallow defects in narrow 
gap semiconductors which can be modeled by the Kane Model* 1 ®*. Much of this 
materials comes from the book by Bastard* 19 *. The essential fact to note is that band 
mixing causes the situation to be much more complication than a hydrogenic model 
- except in the crudest approximate. In any case, for shallow impurity levels, we 
always assume a long range potential with perhaps a central cell connection if some 
refinements are necessary. 

The main purpose of including Appendix C is to give one a feel for the complexity 
of the problem and to show why, since the binding energy can never be appreciable, 
that one often assumes total ionization of the shallow levels. 
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B. Other Complexities^ 20, 21 ) 

We have already mentioned that the donor ionization energy may be so small that 
it is only sensible to assume its total ionization. 

In addition with a large density of states, the donor levels may split into a band 
which may merge with the conduction band. The merging of the donor levels with 
the conduction band seems to be a common occurrence for Mercury Cadmium 
Telluride with x < 0.2. There are other reasons for the merging of the donor levels 
with the conduction band. Donor levels can be screened by the electrons in the 
conduction band, and this screening can decrease the binding energy (Neumark). In 
addition, when the binding energy gets small enough, the donor wavefunctions 
may hybridize with the band functions and again we get merging (Lowney). 


12 



TABLE 1 


Types of Point Defects 


Substitutional Atoms 
Interstitial Atoms (Self or Foreign) 

Vacancy 

Antisite Defect in Compound Semiconductor 
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TABLE 2 


Materials 


MCT 

HgCdTe 

Hgi.xCd.Te 

MZT 

HgZnTe 

H Sl-x Zn x Te 

MZS 

HgZnSe 

H 8l-x Zn x Se 
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TABLE 3 


Brief Properties of Shallow and Deep Levels in II-VI Compounds 



DEEP LEVEL 

SHALLOW LEVEL 

Definition 

Originates from 
short range potential 

Originates from long range potential 

Wave Functions 

Typically mix with 
Valance and Conduction 
band 

May mix only with adjacent band to 
level 

Energies 

Often have ionization 
energies comparable to 
half of gap energy - 
not a Hydrogen like 
spectrum - use tight 
binding approximation 

Usually have ionization energies much 
less than gap energy - have a Hydrogen 
like spectrum as modified by dielectric 
constant effective mass, or other 
complications 

Properties 

Act as recombination 
centers, compensators, 
and facilitators for electron- 
hole generation 

These supply carriers and hence control 
conductivity 

Complications 

Hard to identify 
appropriate models 

Questions of screening (which may be 
site dependent) and merging with band 


Non parabolic bands 

Non parabolic bands and band mixing 



Polar nature of lattice, degeneracy or 
near degeneracy of bands can also cause 
complications 

Examples 

Vacancy, substitutional 
atom, antisite, interstitial 
atom 

Hg vacancy. In on metal site plus many 
others as predicted by periodic table 

Historical Theory 

Slater-Koster 

Kohn-Luttinger Effective Mass 
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TABLE 4 


Summary of Calculations which we Have Done on Deep Levels 


No Relaxation 
MZT 


cation-site-s-like - Te (antisite), C, At, I 
cation-site-p-like F 

anion-site-p-like Cd, Hg (antisite), Zn (antisite). Mg vs. x 

anion-site-s-like there are no deep levels predicted in the band gap 
several x= 0.3 levels - both s- and p-like 
several x = 0.2 levels - p-like 

cation-site-s-like vacancy level is below the valance band edge about 0.4 eV 

MZSl 


anion site p-like 
cation site p-like 
cation sire s-like 
anion-site s-like 


Cd, Zn (antisite). Mg, Be 

Vacancy vs. x 

Vacancy 

there are no deep levels predicted in the band gap 


With Relaxation (when the parameter a = 0, the use of no relaxation is included). 
MZT 


anion-site 

cation-site 

anion-site 

cation-site 

Vacancy 

p-like 

s- and p-like 
p-like 

s- and p-like 
s- and p-like 

x = 0.15 
x = 0.15 

x = 0.3 
x = 0.3 
x = 0.15 

a = - 

0.4, - 0.2, 0.0, 0.2, 0.4 

MCT 





anion-site 

cation-site 

p-like 

s- and p-like 

x = 0.22 
x = 0.22 

a = - 

0.4, - 0.2, 0.0, 0.2, 0.4 

MZS. 





anion-site 

s- and p-like 

x = 0.08 

(a = 

- 0.4, - 0.2, 0.0, 0.2, 0.4) 
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TABLE 5 


Bond Lengths 

parameter a s a p and p-like deep levels with or without relaxation 


System 

d H (A) 

dj (A) 

a s (eV) 

a p (eV) 

deep levels 
(no relax) 

deep levels 
(relaxed) 

MZT : Hg Xe 

2.74 

2.80 

0.07 

-0.04 

C. B. a 

C. B. 

MZT:Zn Te 

2.74 

2.58 

-0.22 

0.13 

0.9 Eg b 

0.4 Eg 

MCT : Cd Te 

2.54 

2.71 

0.20 

-0.11 

0.98 Eg 

C. B. 


a. C. B. means conduction band resonance. 

b. Eg means the energy band gap. 
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vacancy level 

FIGURE 1 

Deep Levels in MZT (cation site, s-like) 
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FIGURE 3 

Effects of Lattice Relaxation on Deep Levels in MZT (anion site, p-like) 
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FIGURES 

Effects of Lattice Relaxation on Deep Levels in MZS (anion site, p-like) 
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B. Green's Function^ 1 ) 


1. Summary of some standard Results^) 

The results we list can be found, together with derivations, in the book by 
Economou. 


The Green's function can be defined by 




=x 


n> <n 

" E-E„ 


(Bl) 


Define 


-* — * lim -* — * 

G±(r,r’;E) = + G (r, r’; E± is) 

s — » 0 


where G (r, E) = (r I G (E) I ?) 

One can then show that the local density of states can be derived from 


P (E) = T -Im 

K 


Tr 


(G ± (E) )] 


(B2) 

(B3) 


(B4) 


where we have used: 


G* (r, r^ ; E) = 
* 

<M r ) <K ( r ~) 

n r _ tt 


lim <Mr) (Mr 1 ) 
s->0 + E ± is - E n 

Ti*^8(E-Ej* n (?)♦„(?) 


(B5) 


Consider an unperturbed system with Hamiltonian H Q and a perturbed system with 
Hamiltonian Hj. Let the corresponding Green's function be G 0 and G where 


G 0 (z) = (z-H 0 )- 1 

G (z) = (z - H)' 1 (B6) 
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From this one can derive the Dyson equation 


G - G 0 + G 0 H i G 


And if one defines T matrix as 


T = H, G Go 1 
one can show 

G(z) = G 0 (z) + G 0 T G 0 , 

thus from the T matrix, we can determine the Green's function. 


(B7) 


(B8) 


(B9) 


We now apply this to a crystal with an unperturbed tight binding Hamiltonian 
given by 


H 0 = e 0 Ii|l> <1 |+I hl l>Vi^<m 


(BIO) 


where we will restrict Vf . m to be ^ for 1 and m being nearest neighbors and 2 
otherwise. 

Note II > is the atomic-like orbital, or Wannier function, centered at lattice site 1 . 
We will consider only one band. The Bloch and Wannier functions are related by; 


1 — ► 

1 1 ^ y"* i * .. 

Ik > = —?= Xre 1 1 > 

Vn 1 


(Bll) 


Since H r 


k > = E (k)| k> , we can show; 




(B12) 


the typical tight binding result. Knowing the exact energy eigenvalues we can then 
evaluate the unperturbed Green's function. 


27 



To show how this is useful we consider a simple model of a substitutional impurity 
in a lattice as given by Economou. In Economou's model the substitutional 
impurity has as a perturbation 


H 1= II >e<l 


From this, one can obtain the T matrix as 


(B13) 


T = 1 1 > 


l-eG 0 (l,l;E) 


<1 I 


(B14) 


where G 0 (1, 1; E) = < 1 I G^ I 1 > 

The Green's function is then given by; 


G = G 0 + G 0 I 1 > 


l-eG 0 (l,l) 


<1 I G r 


(B15) 


Then the poles of G yield the discrete eigenvalues (Ep) of H via 


G 0 (r,r,E p )=V E . 


(B16) 


One might also mention that an ideal vacancy is modeled by assuming e — » <» or 


G 0 (l ,1 ,E ) =0 


(B17) 


which then excludes the side 1 from occupancy. 
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2. Analytical Results 


The Green function G 0 (E - H 0 ) has a singularity at E = E n £• , the eigenvalues of 
H q . The singularity may be by-passed by including complex values 


Go(E)= L 


e-^o E - H 0 ± i e 


But 

L 1 1 

= p — q: i 7r 5 (E) 

e-»o E±i e E T 

where the principal value is defined by 


f h 

r*- 5 f b 

| f W = L 

| f (x) dx ( J f (x) dx 

1 x-x G S^o 

1 X - X 0 1 X - X Q 

J a 

j a *' x o + 5 


Thus 


(B18) 


(B19) 


(B20) 


Go(E) = P Fir :jr 

and the first term on the right hand side can be rewritten (formally) 


± I 5 (E ' - H ) 

G c (E) = P I — _ £ , dE- T(in)8(E-H e ) 


which is another form of (B5). 


(B21) 


(B22) 


29 



Away from the singularities we can write 


G 0 (E) = G 0 (E) = 


8 (E ' -H ) 
E-E' 


dE’ 


Using the identity 


8(E’ -H°) = £lnk>8(E' -E nk )<nk 

nk 


_ , I n k > 8 (E ' -E.)<nk I 

G » (E) ^I — 


(B23) 


(B24) 


(B25) 


Using I 1 > to denote Wannier state localized at site 1. 


„ <1 I n k > 8 (E ' - E.) < n k 1 1’> 

<1 IG C (E) ll'> = X , " k dE' 

n k k - k 


<i IG„(E) H'>=X nk 


<1 I nkxnk 1 1 ’ > 


E-E 


n k 


(B26) 


(B27) 


Finally, it is very useful for relaxation calculations to derive an expression for the 
force in neighboring ions. 

Formally, the total energy can be written as 


/ 


E tot = I Ep(E)dE 


(B28) 


where p (E) is the density of states. 

The force on one of the nearest-neighbor atoms can be represented as 
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f aE tot_ g 

x 3 x 3 x 


| E p (E) 


(B29) 


and x is along the defect-host bond direction. More precisely, this doesn't treat the 
core electrons which give repulsion, so this is only for the attractive electrons, F x • A 
basic result is (Eq. B4) 


P (E) = ~F “ Im TrG (E) 


(B30) 


Thus 


/• E F 

F* = + I E — Im [t r G + (E)l 
x 3 x I it L J 

» - oo 


(B31) 


Dyson's equation works for G(z), so it certainly works for G + (E) which is just a limit 
ofG(z). Thus 


G = G 0 + G 0 VG 


(B32) 


*C-_ G /V G+ 3G 
3x °3x ° 3x 


(B33) 


( i - g ° v )H =g °!t g 


(B34) 


which means 


||.(I.G o V)- 1 G 0 ^G 


(B35) 


From Dyson's equation it also follows that 
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(B36) 


(B37) 


(B38) 
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3. Determination of G°(E) - Vogl Band Structure Theory* 2 ) 

To calculate the defect levels from Eqn. (6), we need to calculate G s / G p first. In the 
Appendix B2, we have shown for E in the band gap as for defect levels (i.e. with E * 
E^ k ), we have 


G C (E) 


r 

8 (E ' -H°) 
E-E’ 

j 


dE 


(B39) 


We have also noted that the spectral 


density operator 8 (E ' - H ) can be written as 


8 (E ’ -H°) = X 1 nk>8(E-E n £)<nk I 
n iT 


(B40) 


Where E nk and I n k > are the eigenvalues and eigenvectors of the host crystal. In 
turn, E n £- and I n iT > can be determined from the model of H 0 . 


The model of H 0 which has been used is based on the model of Vogl. The 
Hamiltonain of the perfect crystal has the form: 


H 0 = ]jr|liaR>Ej J <iaR 1+ licR’>Ej<icR' ij 
+ £ I i a R > Vj j< j cR+d I + h.c. 


(B41) 


where i = s, p x , Py , p z or s* labels the orbitals, a and c denote anion (like Te' ') or 

cation like (Hg ++ , Zn ++ ), and h.c. denotes the Hemition conjugate. The states 
li a R > and li c R' > are localized orbitals centered on the R and the cation at 

R + d , respectively, where d is the nearest neighbor distance. 

The main points of this model as summarized by Vogl are 

(1) The chemistry of sp bonding is retained. 

(2) The diagonal matrix elements of the model are related to the atomic 
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energies of the chemical constituents. This is really what allows one 
to predict chemical trends. 

(3) The off-diagonal matrix elements scale as the inverse of the bond 

length. 

(4) The number of parameters in the model is minimal. 

(5) The theory reproduces both valence and conduction bands. 
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4. Determination of V - Hjalmarson's Deep Level Theory^ 3 ) 

From Eqns. (6), (7) in section HA, we see we need the perfect crystal Green's function 

G° and the defect potential V to determine the energy band. In the representation 
with symmetry appropriate to the symmetry of the point defect, V s and V p are 

matrices of a certain size, arising from Bloch diagonizalition. With the five atom, 
central atom and nearest neighbors, cluster, V s and V p are 5 x 5 matrices. If we want 

to include the 2°^ neighbor in the calculations, we need to consider a 17 atom cluster 
which includes a defect, 4 nearest neighbors, and 12 second neighbors so the 
matrices are 17 x 17. 

In Hjalmarson's deep level theory, four assumptions, described in detail in 
Hjalmarson's dissertation, are given. The assumptions are 

(1) Only consider the central cell potential (ignore the Coulomb potential). 

Thus this theory can only be used to predict the deep levels. 

(2) Use the nearest neighbor approximation so V s and V p are 5x5 

matrices. 

(3) Assume the off-diagonal elements of V are equal to zero and this 

ignore the effects of lattice relaxation. 

Therefore, the form of V is diagonal and can be written: 

V=XiVi = Zi liao>Uj<iao I (B42) 


where the U i are the diagonal elements of V . and the basis states I i a o > are 

orthogonalized symmetrized Lowdin orbitals which are linear combinations of sp 3 
hybrid orbital. This assumption has been called the ideal defect assumption. 


(4) Assume the diagonal elements of the defect potential Vj can be 
calculated with the equations 


U 


S 



(B43) 


u p =P p (e^-eLJ 


(B44) 


where E imp and E host are the atomic orbital energies of defect and host atoms 
repsectively, with (5 S and (5 p being proportional constants. Hjalmarson indicates that 
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(^imp " ^host) changes by a constant factor on going from are free ions to ions in a 
solid, so that this approximation puts in the correct chemical trends. Hjalmarson 
and other discuss empirical ways of fitting the P's. 

Because the V s , Vp are diagonal, it is easy to show that the defect eneriges are 
determined by 

U = 1 /s and 

G aa 

Where G aa and G^ a 
G s ,G p respectively. 


U P =/ G P 

aa 


(B45) 


are the diagonal elements of the Green's function matrices 
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C. Details of Shallow Impurity Levels^ 19 ) 

Shallow Impurity Levels and the Kane Model 

We start by reviewing the Kane k • p three band model. Kane diagonalized the 

Hamiltonian exactly for three bands (conduction, light hole, and split off). We write 
the basic Hamiltonian as 


H = 


2m, 


+ V(r) 


where p is the momentum operator, nr 0 
potential. Using Bloch's theorem 

W(r) = VQ e Ur(r) 


(Cl) 

is the electron mass and V(r) is the crystal 


(C2) 


where Q is the volume of the crystal, k is the Bloch wave vector and Uj^(r) is 

periodic between unit cells. The time independent Schrodinger equation can then 
be written 


H k + 


Tr ir 


(k - k 0 ) - p + (k 2 - k^) 


o m. 


2m, 


U„ k W 


= E„(k)U„ k (r) 


where 


2 2 2 

p 2 - -ft _ _ 

H k =^— + V(r) + — k 0 .p + 2 

° 2m m c 0 v 2m 


The band index is n. 


(C3) 


(C4) 


The equations for the k • p representation are obtained by expanding the U n ^* in 
terms of a set u nko 
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U„lM = I„.K nn .U„. Eo (r) 


(C5) 


K jT /q » U "SW U "So« dv = 5 r> n 


(C6) 


* 

over a unit cell then combining the last two equations, multiplying by ^nk Q 

/ 

and integrating over a unit cell, we find: 


E " ( k ° )+ ^ H 1<2 ' k “ 1 

A All o 


^nn + m ( k " k o) • Pnn / K nn 
o 

= E„(k)K„„, 


(C7) 


where 


Pnn ' Q 0 ho U nk 0 P U n k 0 dv - 

The spin-orbit interaction must also be introduced, 
interactions. This amounts to adding in 

H so = — \~2 ( ax ^v). p 

4m 0 C 


(C8) 

We ignore other relativistic 


(C9) 


where a is the Pauli operator. 
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Including the spin-orbit term one finds: 


2 , 2 


In' - It E n Co) - E„(k) ♦ ] s„„. + A- k . p nn . 

+ — — — < n l(oxVv) . p I n ' > 

m o 4m 0 C 2 

tlk ^ <n l(oxVv) ln'>}K nn . = 0. 


o 4m 0 C 


(CIO) 


Where the last term comes from going from V n £ to U n £ and is very small and 

will be neglected from here on out. We have chosen the T point where k Q = 0 for 
Eqn (10). 

The approximation is often used that U n0 is restricted to states of r 6 , T 7 , and r 8 

symmetry - in the usual sense. These states include four bands (the conduction 
band, the light hole band, the split off band and the heavy hole band). 

The basis functions for the bands, which diagonalize the spin-orbit interaction, will 
now be explained. We start with the functions ls>, lx>, ly>, lz> which transform 

just like atomic s and p functions under the action of the tetrahedral group. 
Including spin, this gives 8 functions: 

Is1>, lxT>, lyt>, lzT> and I sl>, I xi>, I yi> and I zi > . We then form linear 
combinations which diagonalize the spin orbit interaction. 

In order to do this, we form linear combinations for which J = L + S and its 

projections J z are diagonal. We summarize this in a table below, together with the 

band edge energies. This is a problem which is solved by Clebsch-Gordon 
coefficients. The U n below are the periodic parts of the Bloch functions. In what 

follows we follow Bastard (1988) very closely, but add some details for clarity. 
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Linear Combination 


Energy of 
Band Edge (H=0) 
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We will display the matrix of the Hamiltonian later, but for now we give some 
results. The secular equation for the energy is; 


E 3 + (A-E^E 2 -(E g A + P 2 k 2 )E-|AP 2 k 2 = 0, 


(Cll) 


from the three coupled bands. The heavy hole band decouples and comes out 
dispersionless. However, coupling of the heavy hole band to higher bands leads to 


E hh “ 


2 2 
ti k 


2m 


hh 


(Cl 2) 


The momentum matrix element P and the spin-orbit splitting parameter A are 
defined by 


P = < s I PI x > 

m „ x 


(C13a) 


and 


A = -3 i i 


2 2 
4m 0 C 


< X 


I (vv X p) 2 


I z > . 


(C13b) 


Our next job is to work out the matrix elements. We drop the k-dependent spin- 
orbit term as well as the relativistic corrections. 


2 2 

M ij =(E i (o)) + ^^8 :; +— k . Pl 


2m o m o “ 


m o 4m 0 C 2 


i I (oxVv)*P 


J > 


The results are the following: 


(C14) 
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, 1 

2 

3 

4 

3" 

6 

7 

8 

1 

2 uric 

0 


>*• 

PAfe + 

0 

3^1 

-Jink. 

2 

0 

Ml 

2 m* 

■J3^ + 

-JfPfcfez 

0 

fhk- 

MK 

iff ^2 

3 

-3^ 


E 


(9 

0 

0 

£ 

4 

-j|P/ik + 

B 

£ 

-etEfe* 

c * aw. 

7) 

0 

0 

0 

5 

m. 

0 

i 

0 

(9 

Hg 

0 

0 

0 

6 

<9 

PAfe + 

0 

(9 

0 

cM 

L ' 2m, 

0 

* 1 

0 

7 

js^» 

Jink. 

0 

0 

0 

0 

1 

KSM 

0 

8 

-3nK 

i* k z 

0 

0 

0 

0 

0 

-£* - a 

+ 4!if 

•2m, 


























































Shallow Defprts 


We want to describe the behavior of electrons moving in imperfect III-V and II- VI 
crystals with defects. We look for solutions of the time independent Schrodinger 


[ H o + ^(r)]v(r)=e\|/(r) 


(Cl 6) 


where H 0 i S the Hamiltonian of the perfect crystal and <|> (r) is the non-periodic 

potential arising from the defect. For shallow defects, <|> (r) is slowly varying at the 
scale of the unit cell. 


We first prove some initial results. Consider the integral 
I = L f ( r ) P ( r ) dV 


(Cl 7) 

where p (r) has the periodicity of the lattice. 


p(r + R) = p(r) 
and R = X n i e i 


(Cl 8) 

(Cl 9) 


with n 4 ( i 1, 2, 3 ) being integers and e, the fundamental translation vectors of the 
lattice. / is required to fall off quickly enough at large distances for I to converge and 
p (r) is periodic. It can be expanded in a Fourier expansion using only the reciprocal 
lattice vector R . 


p(r) = X£p(K)e iK#r 
where by inversion 


(C20) 
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p(K) = l/ 0 PW«'‘ e ' rdv 

= 75 -/ p(r)e“^* r dv, 

ii o J n 0 


(C21) 


where £2 is the volume of the crystal and £ 2 0 is the volume of the unit cell. 


An approximate form for the integral below can be derived. This form of integral 
occurs later and we need to approximate it. 


I 


= { a f (r)p(r)dv 


= | p(K) f^/(r)e iKtr dv. 


(C22) 


Since the volume of the unit cell can be decomposed into identical stacked unit 
cells, we can write 


J a /(r)e lK ' r dv=£ J qj / (R| + r) e‘ K * r dr , 

Ri 


- i K • r 


i K R 

where we have used e = 1 . 


On the scale of a unit cell, we assume / (r) is slowly varying so 




(C23) 


Thus 


f / (r) p (r) dv = Irj/ (Rj) f 

J £2 Jf 

= a„8KSR,/(Ri) 


i K • r , 

e dv 


(C24) 
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Thus 


I = Z*P(K)0 0 8°I R1 /(Ri) 

= p(o)Q o X Ri / (Ri) 

= P (o) f a f (r) dv 


But P ( °)=jJ n P« dv ,so 


we find 


l£ Mn P(?)dv )(i /(?)dv )- 


We have assumed the Bloch states 1 

n k > are normalized as per 

S "n = A| u’ k f)U„. k dv 

Jn 

(whole crystal) 

"s;/n 0 u -« u " k dv 

(unit cell) 

To solve [H 0 + <J> (r)J y (r) = E \|/ (r) 



(C25) 


(C26) 

(C27), 


we assume <J> (r) couples the T, , r 7 , r R . . 

w v 678 states among themselves. 

a Kane like theory we look for a solution in the form 


Again, following 




(C28) 


where the / x (r) 's are the envelope functions which vary slowly. The U lo 's are the 

, _ 1 i k • r 

T point Bloch functions for the eight bands. When (|) (r) = 0 , we have * l ^Q 6 
The wave function \|/ (r) is normalized so; 
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■L 


2 

V(r) I dv 


m J^ / m /]U mo U lo dv 

we^aup 1 ^ E<5S and ^ f ° r mtegrals of P eriodic and slowly varying functions. 


1 Zi m C 1 C m —jf m f I dvju mo U lo dv 


So we find ; 

i=i 'N 2 Jj/i 


dv 


(C29) 


And again since we want X |c, | - 1 , the envelope function needs to be normalized 
to 1. 


Thus 


| j/,f dv = i 


In summary we have: 

p 2 /_ _ . 

H ° = 2^T +V(r) + — T-2 (ox vvj.p 

° 4 m o C 

V(r)=sf. 1 C 1 /,f)U lo (?) 

[h o + ^ ( r )] V (r) = e \j/ (r) 


(C30) 


(C31) 

(C32) 

(C33) 


and /, is normalized by Eq (30). We now substitute Eq (32) into (33) and use (31). 
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We use 


P 2 (/l^lo} = ( P 2 /l)U| + 2(p/|)(pU|)+ /|(p 2 U,) . (C34 

LetH o be 8 iven by E q (31), and as usual we will neglect the "k - dependent" part of 

thespin orbit interaction. We also set the notation so the spin orbit term is buried 
in J-T , 


H oU lo =e lo U lo 


So we get 


0 = I 1 8 =1 c i{( e io-e)|/m/iU, o uL dv 
+ | Q Urn ° u i°/ m ^ r /i dv 

+ m J n P bljo) • (/mP/i)dv . 

0 •'Ll 


+ l Q fL l K'>/,uLu lo dv 

Using again Eq (25), we can simplify several integrals. We get 

0 = S?, 1 c i | ( e rc)^|/^/|dv|u^ o u, o dv 

2 

+ n/^ m ^^ ldV j U mo U lo dv + ^-j/m < t>(r)/lju mo U lo dv 

+ P ml°|/mP/l dv )=0, 

where 

Ilf • — 

Pml ~ m Q Q J Q Umo P U io dv / 


(C35) 


(C36) 


(C37) 


(C38) 
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and we can use 


I 


u modv = ns m . 


Therefore, we write 

0=Z®,, C l(( e lo- e )|/m/l dvS lm + |/m^-/l dv5 l 
/mOW/,dv8 lm+ P ml .f45/ | dv) 


m 


We can rewrite this as 


°=| si dv /L 


»sr.i *i 


m 


Ej - 8 + 


2m, 




+ Pm 1 • / 1 / 


where due to periodicity, we can rewrite Eq (38) as 
1 


P ml = 


m^J U « PU '- d *’ 


Non vanishing solutions for / exist if /' is an eigensolution of the problem: 


(» + ♦) 




r\ 

/ 2 


h 

• 

= £ 

* 



- 00 


(C39) 


(C40) 


(C41) 


(C42) 
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where and the matrix D + <J) is given on the following page. We have 

dropped the free electron term and continue to neglect relativistic corrections and 
the 'k dependent' part of the spin orbit term. 
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The matrix for D + <{> can be written: 



1 

z 

3 

4 

5 

6 

7 

8 

1 

<#>-£ 

0 

-fm 

Arr 

VP t 

0 

km 

-fi PP- 

2 

0 


Am, 

-Jim 

£ 

pp- 

Jim 

km 

3 

-Jim 

JS r - 


0 

0 

0 

0 

0 

4 

-Am 

-Mr Pi 

0 


0 

0 

0 

& 

5 

TP- 

0 

0 

0 


0 

0 

0 

6 

0 

PPr 

0 

0 

o 


0 

0 

7 

Am 

J%?P- 

0 

0 

0 

0 

0- A 

0 

8 

-Jim 

£v 

J3 r Z 

0 

0 

0 

0 

0 

0-A 

-t.-i. 
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The periodic part of the Bloch functions going along with f x are 


Ui 

U| = i I s T > 
U 2 i \si> 


U 3 = 

-VF izT>+ 

76 ' lx+ 

i y) i > 

u 4 = 

■7e 

|c^ko 

1 

A 

1 z i > 



■iy)T> 



X 

-Kf 

II 

so 

D 

- i y) i > 


IV 

• *— 1 
+ 

X 

-f? 

II 

> + _ r 1 

73 

z T > 

U 8 = 

-w llx - iylT 

> + ~r 1 

73 

z i > 
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Eq (42) can be reduced to an effective 2x2 system starting with the third row, we get 


V 2 P 

3 p Pz/i + ^P-/2 + (-e-eo +< t>)/3 = 0 ^ 


or 


/ 3 “ 


( e o + £ ’ ♦) L 


' VT PPz/l + V3 P - * 2 


In a similar fashion, from the rest of the rows, we can derives 


(C43) 



(C44) 

(C45) 

(C46) 

(C47) 

(C48) 


If we substitute Eqs (43) - (48) into the two equations we get from the first two rows 

of Eq (42), we get our effective 2 x 2 Eq (42). This, in effect, projects the Eq into the r 6 
band. 
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The resulting effective 2x2 equation can be written 




/ 1 
f 2 


= £ 


/ 1 
f 2 


where 



CM 

k 

_1 2 

( l-l 

*d= 

3 Pzle + eo-<l>| 

P z+ p ?+ 

(e+e 0 -<|)) p- 

2 

i r 1 

2p2 i 

\'i 

rP 

> 

-l e + e o-4>J P + 

+ 3 P -( e 

+ e 0 +A-<| > ) p + 


+ yPz( e + e 0 + A -<1>) P z +<t> 


(C49) 


(C50) 


a/2 2 

h -=t p 


. {(e+e 0 + A-(|)) -[e+e 0 + A-<|>]^p z -p z |[e+e 0 -<|>] -[e+ e 0 + A - <t>] j 


l P - 

(C51) 


Hvdrogenic Donor Problem 



K = static dielectric constant of semiconductor. 


(C52) 


For almost all r, (<[> (r)| < < e o 
Expand to lowest order in 


-<1> + e 


-<)>+£ 
e 0 + A 
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We obtain 


r 6 


(C53) 


„ „o 2P 2 2 p 2 e 2 2P 2 x - 

H d =H d 2" e P TP + — Pz^WPz 

3e 0 3 (e 0 + A) 3e 0 

P 2 p 2 

+ -yP + 4>(r)p.+ — p.<j)(r)p + 

e o 3e 0 

P 2 2P 2 

+ - Pz ♦ O') Pz + -p. <t> (r) p + 


3(e 0 + a) 


3 (e Q + a) 


4l 2 

H nd = ^P 


( £ o +A ) 2 


P-^WPz'Pz^W P- 


(C54) 


With the zero order term only get Hydrogenic levels 


m 


a„ = 0.53 A x K x — ° 


m 


m, 


R 0 =13.6eV x —L 
K m. 


(C55) 
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